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We propose a method that allows for the efficient solution of the three-body Faddeev equations in 
the presence of infinitely rising confinement interactions. Such a method is useful in calculations 
of nonrelativistic and especially semirelativistic constituent quark models. The convergence of the 
partial wave series is accelerated and possible spurious contributions in the Faddeev components are 
avoided. We demonstrate how the method works with the example of the Goldstone-boson-exchange 
chiral quark model for baryons. 
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I. INTRODUCTION 

The Faddeev approach has for a long time been a very 
efficient tool for solving three-body problems in nuclear 
physics. Especially in the case of short-range forces the 
Faddeev equations can nowadays be solved to a high 
precision both for bound and scattering problems. The 
Faddeev theory offers a number of advantages over other 
approaches, such as solving the Schrodinger differential 
equation. In particular one can immediately include all 
relevant channels with the proper boundary conditions; 
in the cases of either two or three identical particles, a 
straightforward scheme for antisymmctrization is imple- 
mented; and finally a wide class of (local and nonlocal) 
two- and three-body interactions can be dealt with. 

However, with the Faddeev integral equations one en- 
counters difficulties whenever the interactions are long- 
ranged. There is the notorious problem of solving the 
three-body Coulomb problem. Similarly, the solution of 
the constituent quark models poses practical difficulties. 
Irrespective of the type of confinement, one is confronted 
with complications due to the fact that the potentials are 
rising infinitely, as a function of inter-quark distance. 

In the past the Faddeev theory has already been suc- 
cessfully applied to the calculation of nonrelativistic con- 
stituent quark models. In this case, when the confine- 
ment in general remains relatively weak, a direct solution 
of the usual Faddeev equations can be achieved with a 
sufficient accuracy, and a variety of results for baryon 
spectra has been obtained [Q-Q. 

The situation becomes much more difficult for semirel- 
ativistic constituent quark models (i.e. with the kinetic- 
energy operator taken in relativistic form); in this case 
the confinement is considerably enhanced, to a strength 
practically consistent with the string tension of quan- 
tum chromodynamics. The convergence of the partial 
wave series gets slow, a relatively large number of chan- 
nels needs to be included, and one may also easily pick 
up spurious contributions in the Faddeev components, 
which get wiped out only in the total three-body wave 
functions. The corresponding problems are, of course, 



already there in a nonrelativistic framework and they 
were thoroughly studied (e.g. in Ref. (|]) for the case 
of the harmonic-oscillator potential, which sometimes is 
also used as a simplified model for confinement. 

In the attempt to make the Faddeev approach directly 
amenable to any three-body problem with confining in- 
teractions, be it in a nonrelativistic or semirelativistic 
framework, one of us has previously suggested a differ- 
ent set-up of the Faddeev formalism ||. The proposed 
method exploits the asymptotic filtering property of the 
Faddeev scheme by adopting a different splitting of the 
full three-body Hamiltonian. Already in the nonrelativis- 
tic case considered in Ref. |J this has led to a more effi- 
cient solution of the (modified) Faddeev equations. Here 
we further elaborate on this method and demonstrate its 
efficiency even in the case of semirelativistic constituent 
quark models. 

Usually the full Hamiltonian H of the system under 
consideration is split into a free Hamiltonian (kinetic- 
energy operator) H° and an interaction part H 1 ; the 
latter contains all two and three-body forces. However, 
this is not absolutely necessary. One may equally well 
adopt a different splitting of H , where part of the inter- 
actions are put into a modified Hamiltonian H° and the 
residual ones are retained in H 1 . This turns out to be es- 
pecially useful whenever there are interactions that cause 
difficulties in the original Faddeev scheme. It is usually 
the case for interactions that do not fall off fast enough 
at large distances, e.g., the Coulomb forces and notably 
also infinite confinement potentials. 

In the following section we outline the method that can 
treat confined three-body problems in an efficient and ac- 
curate manner. Beyond solving nonrelativistic problems 
we demonstrate how the method can also be applied to 
semirelativistic constituent quark models, which rely on 
a kinetic-energy operator in relativistic form. In Section 
3 we show how to solve the modified Faddeev equations, 
and we prove the good performance of the method in 
Section 4, taking the specific examples of the Goldstone- 
boson-exchange (GBE) constituent quark model in the 
nonrelativistic and semirelativistic versions. 
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II. TREATMENT OF CONFINEMENT IN THE 
FADDEEV APPROACH 

We start from the total Hamiltonian of a nonrelativis- 
tic or a semirelativistic three-quark system, which can be 
written as 



vp 



(1) 



where H° is the three-body kinetic-energy operator and 
v s = v% + vJ represents the mutual quark-quark interac- 
tions containing both the confinement (vf) and hyperfine 
(Vgf) potentials in the subsystems 5 = a, (3, 7. For the 
moment we leave out genuine three-quark forces, which, 
however, can easily be included into the formalism and 
the method we follow. In the nonrelativistic case we 
may express the kinetic-energy operator by four equiv- 
alent forms 
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i.e. either through individual particle momenta ki in the 
center-of-mass system or in terms of relative momenta 
ps and qs conjugate to the usual Jacobi coordinates xg 
and y$, respectively (5 = a,/3, 7). In Eq. (j^), m, denotes 
the individual particle mass, fi$ the reduced mass in the 
two-body subsystem 5, and Ms the reduced mass of this 
subsystem with the third particle S. In the semirelativis- 
tic case the kinetic-energy operator takes the form 
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where again ki are the individual particle three-momenta 
in the frame with total three-momentum P = 5^=1 ^ = 
0. We note that a Hamiltonian as in Eq. (uh together 
with the relativistic kinetic-energy operator ([|) repre- 
sents an allowed mass operator in the point-form formal- 
ism of Poincare- invariant quantum mechanics, irrespec- 
tive of the dynamical origin of the interactions pQ ]. 

In the conventional Faddeev treatment the total 
Hamiltonian (jl]) is first split into the free Hamiltonian 
H° and the interaction Hamiltonian H 1 = v a + vp + v 1 . 
Then a decomposition of the total wave function 
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is carried out such that each of the three components 
satisfies the Faddeev integral equations of the type 



|VO = G a (E)v a (\iPp) + |V 7 )) 



(5) 



with a, P, 7 a cyclic permutation. Herein the so-called 
channel resolvent 



G a (E) = (E-H <) 



(6) 



occurs, which contains the interactions in subsystem a 
but is otherwise characterized by a free motion of the 
third particle relative to the two-body subsystem. This 
scheme works well for short-range interactions and makes 
it possible to easily incorporate the correct asymptotic 
behavior. Strictly speaking the standard Faddeev scheme 
applies only for potentials falling of fast enough at large 
distances. When the interactions are long-ranged, how- 
ever, there is never a free motion of the third particle 
relative to the two-body subsystem. This makes it nec- 
essary to modify the Faddeev formalism. Otherwise one 
risks unpleasant properties in the Faddeev components. 
In particular, for infinitely rising potentials spurious con- 
tributions are picked up and also the partial-wave series 
becomes slowly convergent. 

One can circumvent these difficulties by modifying the 
decomposition (Q) such that all the long-range poten- 
tials are included in a modified channel Green's operator. 
Specifically, in our case at least the long-range parts of 
the confinement interactions in all subsystems a, f3, and 7 
should be included in the modified channel resolvent M . 
One can attain this goal by adopting a different splitting 
of the total Hamiltonian into 



where 



H = H c + v a + vp + v n 



(7) 



(8) 



contains, besides the kinetic energy, the long-range parts 
tig of the confining interactions Vg in all subsystems. The 
potentials vg are the residual interactions containing the 
hyperfine potentials and the short-range parts of the con- 
finement. 

Based on Eqs. (0) and (||) we now decompose the total 
wave function into 

l*> = IVU + IV73} + |V7>, (9) 
where the modified Faddeev components are defined as 

|VO =G c (E)i a \t>) (10) 

with 

G C {E) = (E-H C )- 1 . (11) 
They fulfill the integral equations 

$ a ) = G c ct {E)v a {\j>p) + \i> 1 )), (12) 

with a, 0, 7 again a cyclic permutation. As compared to 
Eq. (0) the new channel resolvent gets modified to 



G° a {E) = {E-Ht-Vo)- 1 



(13) 



It exhibits just the desired property of including the long- 
range confining interactions in all subsystems a, (3, 7. 
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Only the short-range potential v a remains in the modified 
Faddeev equations (|l2|). Specifically, since now G c a con- 
tains also the long-range parts Vp + of the confinement 
interactions in channels (3 and 7, the dependence of the 
component \ip a ) on the Jacobi coordinate ys can never 
become a free motion. Rather the proper confinement- 
type asymptotic conditions are imposed on \ip a )- As a 
result, spurious contributions are avoided in the indi- 
vidual Faddeev components, and at the same time the 
partial-wave expansion converges much faster. 

The splitting of the interactions in Eqs. ([?]) and (g) has 
to be done with care. In general, the interaction parts 
put into H c must not produce any bound state. Other- 
wise the proper behavior of the Faddeev components \tjja) 
would again be spoiled. Suppose the potentials contained 
in H c would produce bound states. Then at the cor- 
responding energies, the resolvent G C (E) would become 
singular. Consequently, according to Eq. (|To|), any large 
Faddeev component \tpa) could be generated even if the 
full solution |W) remains infinitesimally small. Therefore, 
besides the true physical solutions of the Hamiltonian H , 
Eqs. ( [l2l) would also produce spurious solutions associ- 
ated with the discrete eigenstates of the Hamiltonian H c 
jtjj . These spurious solutions would occur for any v a , 
thus having no bearing for the physical spectrum of H. 
Of course, when adding up the three individual Faddeev 
components these spurious solutions would cancel out. 
However, they would cause numerical instabilities in the 
practical calculations. Therefore they should be avoided 
by not allowing H c to produce any bound states. 

In the case of confinement interactions the above re- 
quirement cannot strictly be met, since even the longest- 
range parts of the infinitely rising potential generate 
bound states. However, there is a practical way out: one 
needs to eliminate the bound states generated by H c only 
in the region of physical interest. Outside that domain, 
i.e. reasonably far above the physical spectrum, they do 
not matter. In practice, upon splitting the interactions 
in the Hamiltonian (Q) an auxiliary short-range potential 
is introduced with no effect on the physically interesting 
states. It only serves the purpose of cutting off the con- 
finement interaction at short and intermediate distances 
thus avoiding low- lying bound states of H c . In the solu- 
tion of the Faddeev equations below we take a Gaussian 
form for the auxiliary potential as used already in Ref. 



III. SOLUTION OF THE MODIFIED FADDEEV 
EQUATIONS 



We solve Eqs. (|12|) along the Coulomb- Sturmian (CS) 
separable expansion approach jl^. For a two-body sys- 
tem in an angular momentum state I one makes use of 
the CS functions O] defined by 
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in configuration and momentum spaces, respectively. 
Here, L are the Laguerre and G the Gegenbauer poly- 
nomials, and b is a parameter. The CS functions form a 
complete set 
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where, in configuration space, \n) is defined by (r\n) = 
(r\n) jr. 

In three-body Hilbert space one extends the basis by 
defining the direct product 

\nv) a = {\n) a ®\v) a }, {n,v = 0,1,2,...), (17) 

where the states \n) a and \v) a are associated with the 
Jacobi coordinates x a and y a , respectively. The curly 
brackets stand for the angular-momentum coupling in- 
cluding orbital angular momentum, spin and isospin. 
The completeness relation now takes the form 
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It should be noted that we can introduce three analogous 
bases, which belong to different fragmentations a, /3, and 
7- 

In the practical solution of the modified Faddeev equa- 
tions (|lj) we terminate the sums in Eq. ( |l7|) at some fi- 
nite N. Of course, the basis can always be made large 
enough to produce any desired accuracy. Thereby we 
provide for a separable expansion of the short-range po- 
tential v a in the three-body Hilbert space 
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where v^a — a{nv\v a \n'v') p. Inserting expression (|l9| ) 
into Eq. Q yields 

|VU = G%(E)[l%v a l N $ ) + l%v a ll\^)]. (20) 

By applying the CS states a {nv\ from the left one 
obtains a linear system of homogeneous equations for 
the coefficients of the Faddeev components i> = 

J— a nv 

a {nv\^ a )- 
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A unique solution exists if and only if 
detilG'iE)}- 1 -v} = 0. 
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This condition determines the energy eigenvalues. 

The matrices G C (E) and v have block structure and 
their matrix elements are given by 
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Notice that the matrix elements of the resolvent are 
needed only between states of one and the same parti- 
tion a whereas the matrix elements of the potentials are 
always to be taken between states of different partitions 
a and (3. While the latter can easily be evaluated (nu- 
merically) either in configuration or in momentum space, 
the matrix elements G^ a are more involved. In Eq. ( |2l| ) 
the inverse of the matrix G^(E) is needed and thus we 
make use of the approximation 
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i.e. the inverse of the matrix containing the matrix ele- 
ments of an operator is replaced by the matrix contain- 
ing the matrix elements of the inverse operator. Since 
H c contains all the confining interactions, this approxi- 
mation does not modify the character of the spectrum of 
and, like the approximation in Eq. (|l9|), it becomes 
exact in the limit N — > oo. The matrix elements on 
the r.h.s. of Eq. ( |25| ) can be calculated in a straightfor- 
ward manner, as the matrix element of E and the ones 
of the nonrelativistic kinetic-energy operator are known 
in closed analytical form. The potential matrix elements 
can easily be calculated numerically. In the case of the 
relativistic kinetic energy (^) the matrix elements of the 
corresponding operator are most conveniently calculated 
numerically in momentum space. 



IV. PERFORMANCE OF THE METHOD 

In order to demonstrate the efficiency of our method we 
present nonrelativistic and semirelativistic three-quark 
calculations for light baryons. For the quark-quark dy- 
namics we choose the Goldstone-boson-exchange chiral 
quark model The total quark-quark interaction is 

thus given by 

v = v c + v x (26) 
where the confinement is of the linear form 



The chiral interaction is derived from Goldstone-boson 
exchange and is represented by the spin-spin component 
of pseudoscalar meson exchange (for the corresponding 
formulae see Refs. The parameters of the nonrel- 

ativistic GBE chiral quark model can be found in Ref . M , 
whereas the semirelativistic version is specified in Refs. 
H[§. For the nonrelativistic version of this model the 
solutions have previously been obtained in good agree- 
ment by both the standard Faddeev approach and the 
stochastic variational method (SVM) 14), whereas for 
the semirelativistic model only SVM solutions have ex- 
isted up to now. Here we are completing the solutions 
with the Faddeev integral-equation results also for the 
semirelativistic version of the GBE quark model. 

In the splitting of the short- and long-range potential 
parts in the confinement interaction we introduce an aux- 
iliary potential in Gaussian form such that 



and 



v c (r) = Vo + Cr + a cxp[-(r/r ) 2 
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(29) 



As explained in the previous section, this guarantees that 
the bound states of H c are shifted far beyond the energy 
range of physical interest (i.e. the light-baryon ground 
states and resonances). The parameters of the auxiliary 
potential have been taken as ao = 3 fm~ and ?*o = 1 fm. 
By this choice of the parameter values any bound states 
of H c are avoided below 2 GeV. The values of ao and 
ro also influence the rate of convergence but not the final 
(converged) results. 

With the quark-quark interactions prepared in this 
way we can go ahead and calculate the N and A spectra 
by solving Eq. (|2l]). The results for the ground states 
and the excitation spectra are shown in Tables I- VI for 
both the nonrelativistic and semirelativistic versions of 
the GBE quark model. We have also displayed the rate 
of convergence with respect to the number of channels 
included. The convergence rate turns out to be essen- 
tially the same for the nonrelativistic and semirelativistic 
cases. It is practically not influenced by the strength of 
confinement, which is more than three times stronger in 
the latter case. One observes that only a few angular- 
momentum states are needed to attain satisfactory con- 
vergence. All this is essentially due to the modified split- 
ting of the Hamiltonian according to Eq. (0). The con- 
vergence is faster for the A ground state (Table V). Also 
for its first positive- and negative-parity excitations only 
a rather small number of channels need to be included. In 
case of the nucleon and its resonances at most 10 channels 
are necessary. For the separable expansion of the poten- 
tials, according to Eq. (flS|), about 20 terms are perfectly 
sufficient. 

In all cases satisfactory agreement with the results ob- 
tained with the SVM is reached. In the present work 
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we have also succeeded in treating semirelativistic con- 
stituent quark models, whose kinetic-energy operator is 
of the form (||) . Technically the calculation of the semirel- 
ativistic case is more involved, as the matrix elements of 
H° must also be calculated numerically. Otherwise, how- 
ever, the procedure is the same as in the nonrelativistic 
case. 

Finally we remark that we can also calculate wave func- 
tions by putting together the solutions for the Faddeev 
components. The total wave functions are then expressed 
in terms of CS functions. This fact may facilitate their 
use in further applications. 

V. SUMMARY 

We have proposed an efficient method for solving the 
three-body Faddeev equations for the case of infinitely 
rising confinement interactions. The method basically 
consists in splitting of the total Hamiltonian in a way 
different from the usual approach: the long-range parts 
of the interactions are incorporated into a modified re- 
solvent, which takes the place of the usual free resol- 
vent. Consequently the Faddeev decomposition of the 
total three-body wave function is carried out in a differ- 
ent manner, and one derives a set of modified Faddeev 
equations. In solving them one can avoid the difficulties 
associated with long-range confining interactions. 

We have demonstrated the performance of our method 
in the case of the Goldstone-boson-exchange chiral quark 
model both for its nonrelativistic and semirelativistic ver- 
sions. We have obtained a rapid convergence with the 
number of angular-momentum states included, irrespec- 
tive of the strength of confinement. Our results agree 
perfectly with the ones obtained before In par- 

ticular, for the semirelativistic case we have now con- 
firmed within the Faddeev approach the results that have 
heretofore only been available from a variational method. 
It is comfortable to see that the two solution methods, 
which are quite distinct in nature, produce exactly the 
same answers for the light baryon spectra. 
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TABLE I. Masses of the nucleon ground state and the first 
two positive-parity excitations for the semirelativistic (sr) and 
the nonrelativistic (nr) versions of the pseudoscalar GBE chi- 
ral quark model. The convergence with respect to including 
an increasing number of angular-momentum states l x ,l y is 
demonstrated; s and t denote the subsystem spin and isospin. 
In all cases, N = 20 separable terms have been included into 
the expansion. A comparison to the results obtained with the 
stochastic variational method is given. 
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TABLE II. Same as Table I but for the negative-parity 
excitations iV*(1535) — iV*(1520); these two states are degen- 
erate in the pseudoscalar GBE quark model. 



This work was supported by the Austrian-Hungarian 
Scientific- Technical Cooperation within project A- 
14/1998 and by Hungarian OTKA grants No. T026233 
and T029003. 



5 



Channels 


AT(1650) - V*(1700) - V*(1675) 


# 


7 7 


t 


sr 


nr 


i 


Oil 


1 


1963 


i t An 

1749 


o 

I 


10 1 





1649 


1656 




19 1 
1 Z 1 


n 
u 


104 / 


1DOO 


4 


2 1 1 


1 


1 R45 




5 


2 3 1 


1 


1645 


1655 


6 


3 2 1 





1644 


1655 


7 


3 4 1 





1644 


1655 


8 


4 3 1 


1 


1643 


1655 


SVM 


1647 


1652 



TABLE III. Same as Table I but for the negative-parity ex- 
citations N* (1650) -N* (1700) -iV* (1675); these three states 
are degenerate in the pseudoscalar GBE quark model. 
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TABLE VI. Same as Table V but for the negative-parity 
excitations A*(1620) — A* (1700); these two states are degen- 
erate in the pseudoscalar GBE quark model. 
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TABLE IV. Same as Table I but for the negative-parity 
excitations iV*(1680) — A^*(1720); these two states are degen- 
erate in the pseudoscalar GBE quark model. 



Channels 


A(1232) 


A*(1600) 


# 


l X ly S 


t 


sr 


nr 


sr 


nr 


1 


1 


1 


1244 


1249 


1721 


1590 


2 


2 2 1 


1 


1239 


1238 


1718 


1588 


3 


4 4 1 


1 


1239 


1237 


1718 


1588 


4 


6 6 1 


1 


1239 


1237 


1718 


1588 


5 


8 8 1 


1 


1239 


1237 


1718 


1588 


SVM 


1240 


1233 


1718 


1586 



TABLE V. Same as Table I but for the A ground state and 
its first positive-parity excitation. 
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